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■ ABSTRACT 

00 

We perform an analysis of the diffuse low-frequency Galactic components in the Southern part of the Gould Belt system (130° < / < 230° and 
1 ' -50° < b < -10°). Strong UV flux coming from the Gould Belt super-association is responsible for bright diffuse foregrounds that we observe 
, from our position inside the system and that can help us improve our knowledge of the Galactic emission. Free-free emission and anomalous 
, microwave emission (AME) are the dominant components at low frequencies (v < 40 GHz), while synchrotron emission is very smooth and 
1 faint. We separate diffuse free-free emission and AME from synchrotron emission and thermal dust emission by using Planck data, complemented 
• • by ancillary data, using the "Correlated Component Analysis" (CCA) component separation method and we compare with the results of cross- 
^ correlation of foreground templates with the frequency maps. We estimate the electron temperature from Ha and free-free emission using two 
methods (temperature-temperature plot and cross-correlation) and we consistently obtain ranging from 7000 to 2000 K for a dust absorption 
rN fraction of /d = 0-0.5. We estimate the frequency spectrum of the diffuse AME and we recover a peak frequency (in flux density units) of 
25.5 + 1.5 GHz. We verify the reliability of this result with realistic simulations that include the presence of biases in the spectral model for the 
AME and in the free-free template. By combining physical models for vibrational and rotational dust emission and adding the constraints from the 
thermal dust spectrum from Planck and IRAS we are able to get a good description of the frequency spectrum of the AME for plausible values of 
the local density and radiation field. 

Key words. Galaxy: general - radio continuum: ISM - radiation mechanisms: general 
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1. Introduction 

The wide frequency coverage of the Planck^ data gives a 
unique opportunity to study the main four Galactic foregrounds, 
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namely free-free emission, synchrotron emission, anomalous 
microwave emission (AME) and thermal (vibrational) dust emis- 
sion. The different frequency spectra of the components and 
their different spatial morphologies provide a means for sep- 
arating the emission components. In this paper we apply the 
Correlated Component Analysis method (CCA, Bonaldi et al. 
2006, Ricciardi et al. 2010), which uses the spatial morphology 
of the components to perform the separation. The local Gould 
Belt system of current star formation is chosen as a particu- 
larly interesting area in which to make an accurate separation 
of the four foregrounds because of the different morphologies of 
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the components. Gould (1879) first noted this concentration of 
prominent OB associations inclined at 20° to the Galactic plane. 
It was next identified as an Hi feature (Davies 1960, Lindblad 
1967). Along with velocity data from H i and CO combined with 
stellar distances from Hipparchos the total system appears to be 
a slowly expanding and rotating ring of gas and dust surrounding 
a system of OB stars within 500 pc of the Sun (Lindblad et al. 
1997). A recent modelling of the Gould Belt system by Perrot & 
Grenier (2003) gives semi-axes of 373 x 233 pc inclined at 17° 
with an ascending node at / = 296° and a centre 104pc distant 
from us lying at / = 180°. The Gould Belt thickness is 60 pc. 
The stars defining the system have ages less than 30 x 10^ yr. 

The free-free emission from ionized hydrogen is well- 
understood (Dickinson et al. 2003). Yia is a good indicator of the 
emission measure in regions of low dust absorption. Elsewhere 
a correction has to be applied, which depends on where the ab- 
sorbing dust lies relative to the Yia emission. The conversion 
of an emission measure value to a radio brightness tempera- 
ture at a given frequency requires a knowledge of the electron 
temperature. Alternatively, an electron temperature can be de- 
rived by assuming a value for the dust absorption. Values for 
the electron temperature of 4000-8000 K are found in similar 
studies (Banday et al. 2003, Davies et al. 2006, Ghosh et al. 
2012). Radio recombination line observations on the Galactic 
plane (Alves et al. 2012) give values that agree with those of in- 
dividual H II regions, having temperatures that rise with increas- 
ing distance from the Galactic centre; the value at the solar dis- 
tance where the current study applies is 7000-8000 K. 

The spectrum of synchrotron emission reflects the spectrum 
of the cosmic-ray electrons trapped in the Galactic magnetic 
field. At frequencies below a few GHz the brightness temper- 
ature spectral index, y^s, is ranging from -2.5 to -2.7 (Broadbent 
et al. 1989). Between 1.0 GHz and WMAP and Planck frequen- 
cies, the spectral index steepens to values from -2.9 to -3.1 
(Banday et al. 2003, Davies et al. 2006, Kogut et al. 201 1). 

Thermal dust dominates the Galactic emission at Planck fre- 
quencies above 100 GHz. The spectrum is well-defined here with 
temperature Td ~ 18 K and spectral index ranging from 1.5 to 
1.8 (Planck Collaboration XIX 2011). In the frequency range 
60-143 GHz the dust emission overlaps that of the free-free 
emission and AME, making it a critical range for component 
separation. 

The AME component is highly correlated with the far infra- 
red dust emission (Kogut 1996, Leitch et al. 1997, Banday et al. 
2003, Lagache 2003, de Oliveira-Costa et al. 2004, Finkbeiner 
et al. 2004, Davies et al. 2006, Dobler & Finkbeiner 2008a, 
Miville-Deschenes et al. 2008, Ysard et al. 2010, Gold et al. 
2011, Planck CoUaboration XX 2011) and is believed to be the 
result of electric dipole radiation from small spinning dust grains 
(Erickson 1957, Draine & Lazarian 1998) in a range of envi- 
ronments (Ali-Haimoud et al. 2009, Ysard & Verstraete 2010). 
AME is seen in individual dust clouds associated with molecu- 
lar clouds, photo-dissociation regions, reflection nebulae and H ii 
regions (Planck Collaboration XX 201 1). In the present study we 
will be examining the AME spectrum in more extended regions. 



2. Definition of the region of interest and aim of the 
worl< 

The projection of the Gould Belt disc on the sky is a strip that is 
superimposed on the Galactic plane, except towards the Galactic 
centre (Northern Gould Belt) and around / = 180° (Southern 
Gould Belt). In this work we consider the Southern Gould Belt, 




Fig. 1. Orthographic projection (looking towards the Galactic centre in 
the left panel and the Galactic anti-centre in the right panel) of the 
Planck CMB- subtracted 30 GHz channel showing the Gould Belt and 
the region of interest for this paper (defined by 130° < / < 230° and 
-50° <b< -10°). 



which can be approximately defined by Galactic coordinates 
130° < / < 230° and -50° < b < -10° (see Fig. 1). This 
choice gives us a cleaner view of the Gould Belt, because the 
background emission from the Galactic Plane is weaker here 
than towards the Galactic centre. Notable structures within the 
region are the Orion complex, Barnard's arc and the Taurus, 
Eridanus, and Perseus star-forming complexes. All these emit- 
ting regions, including the diflTuse emission from the Eridanus 
shell at -50° < b < -30°, are at a distance of approximately 
500 pc from us and thus they belong to the local inter-stellar 
medium (ISM) associated with the Gould Belt (e.g. Reynolds 
& Ogden 1979, Boumis et al. 2001). 

In Fig. 2 we show the CMB -subtracted Planck data at 1° res- 
olution, compared with the Haslam et al. (1982) 408 MHz map, 
which mostly traces the synchrotron component, the Dickinson 
et al. (2003) Hof map, tracing free-free emission, and the 100 [im 
map from Schlegel et al. (1998), tracing the dust emission. 
The visual inspection reveals dust-correlated features at low fre- 
quency, which could be attributed to AME. There is also promi- 
nent free-free emission, especially strong in the Barnard's arc 
region (towards / = 207°, Z? = -18°). The synchrotron compo- 
nent appears to be sub-dominant with respect to the free-free 
emission and the AME. 

This work aims to separate and study the diflTuse low- 
frequency foregrounds, in particular AME and free-free emis- 
sion, in the region of interest. This requires the estimation of 
the spectral behaviour of the AME (carried out in Sect 4). We 
compare this spectrum with predictions for spinning dust emis- 
sion, one of the mechanisms that is most often invoked to explain 
AME (Sect. 7). Having a reconstruction of the free-free emis- 
sion, we estimate the free-free electron temperature (Sect. 6), 
which relates free-free brightness to emission measure, and in- 
vestigate the dependence of this result on the dust absorption 
fraction. 



3. Description of the analysis 

3.1. Input data 

Planck (Tauber et al. 2010; Planck Collaboration I 2011) is the 
third generation space mission to measure the anisotropy of the 
cosmic microwave background (CMB). It observes the sky in 
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Fig. 2. Gnomonic projections of the region of interest. Top panels: Haslam et al. (1982) 408 MHz map (left); Ha map from Dickinson et al. (2003) 
(middle); and lOOyt/m map from Schlegel et al. (1998) (right) at native resolution. Middle panels (left to right): Planck CMB- subtracted 30, 44, 
and 70 GHz maps at 1° resolution. Bottom panels (left to right): Planck CMB -subtracted 143, 353, and 857 GHz at 1° resolution. 



nine frequency bands covering 30-857 GHz with high sensitiv- 
ity and angular resolution from 3r to 5'. The Low Frequency 
Instrument (LFI; Mandolesi et al. 2010; Bersanelli et al. 2010; 
Mennella et al. 201 1) covers the 30, 44, and 70 GHz bands with 
amplifiers cooled to 20 K. The High Frequency Instrument (HFI; 
Lamarre et al. 2010; Planck HFI Core Team 2011a) covers the 
100, 143, 217, 353, 545, and 857 GHz bands with bolometers 
cooled to 0.1 K. Polarization is measured in all but the highest 
two bands (Leahy et al. 2010; Rosset et al. 2010). A combina- 



tion of radiative cooling and three mechanical coolers produces 
the temperatures needed for the detectors and optics (Planck 
Collaboration II 2011). Two data processing centres (DPCs) 
check and calibrate the data and make maps of the sky (Planck 
HFI Core Team 2011b; Zacchei et al. 2011). Planck's sensitiv- 
ity, angular resolution, and frequency coverage make it a pow- 
erful instrument for galactic and extragalactic astrophysics as 
well as cosmology. Early astrophysics results are given in Planck 
Collaboration VIII-XXVI 2011, based on data taken between 
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Table 1. Summary of Planck data. 



Central frequency 


Instrument 


Resolution 


[GHz] 




[arcmin] 


28.5 


Planck LFI 


32.'65 


44.1 


Planck LFI 


27.'92 


70.3 


Planck LFI 


13.'01 


100 


Planck Wl 


9.'88 


143 


Planck HFI 


7.'18 


217 


Planck HFI 


4.'87 


353 


Planck HFI 


4.'65 


545 


Planck HFI 


4.72 


857 


Planck HFI 


4.39 



13 August 2009 and 7 June 2010. Intermediate astrophysics re- 
sults are now being presented in a series of papers based on data 
taken between 13 August 2009 and 27 November 2010. 

The Planck data used throughout this paper are an internal 
data set known as DX7, whose properties are described in ap- 
pendices to the LFI and HFI data processing papers (Planck 
Collaboration II 2013; Planck Collaboration VI 2013). However, 
we have tested the analysis to the extent that the results will not 
change if carried out on the maps which will be released to the 
public in March 2013. 

The specifications of the Planck maps are reported in 
Table 1. The dataset used for the analysis consists of full res- 
olution frequency maps and the corresponding noise informa- 
tion. We will indicate whenever the CMB -removed version of 
this dataset has been used for display purposes. 

When analysing the results we apply a point source 
mask based on blind detection of sources above 5 cr in each 
Planck map, as described in Zacchei et al. (2011) and Planck 
HFI Core Team (201 lb). Ancillary data have been used through- 
out the paper for component separation purposes, to simulate the 
sky and data, or to analyse our results. The full list of ancillary 
data is reported in Table 2 with the main specifications. 

3.2. Components 

The main diff'use components present in the data are CMB and 
Galactic synchrotron emission, free-free emission, thermal dust 
emission, and anomalous microwave emission (AME). The fre- 
quency spectrum of the CMB component is well-known: it is ac- 
curately described by a black-body having temperature Tcmb = 
2.7255 i^(Fixsen 2009). 

Thermal dust emission dominates at high frequencies. Its 
spectral behaviour is a superposition of "grey-body" components 
identified by temperature T^ust and emissivity index jSa: 

rRj,dnst(v) cx y^^^V[exp(^y//:rdust) - U, (1) 

where k is the Boltzmann constant and h is the Planck constant. 
In the approximation of a single component, over most of the 
sky we have Tdust ~ 18 K and fid of 1.5-1.8 (Finkbeiner et al. 
1999, Planck Collaboration XIX 2011, Planck Collaboration 
XXV 2011). 

The frequency spectrum of the free-free component is often 
described by a power-law with spectral index -2.14 in RJ units. 
A more accurate description (see, e.g. Planck Collaboration XX 
2011) is given by 

rRj,fF(y) cxG(y)x(y/ 10)-^ (2) 



where G = 3.96(T4f -^\y/40)-^-^'^ is the Gaunt factor, which is 
responsible for the departure from a pure power-law behaviour. 
T4 is the electron temperature Tq in units of 10"^ K (re can 
range over 2000-20000 K, but for most of the ISM it is 4,000- 
15,000 K). 

The spectral behaviour of synchrotron radiation can be de- 
scribed to first order by a power-law model with spectral index 
y^s that typically assumes values from -2.5 to -3.2, depending 
on the position in the sky. Steepening of the synchrotron spec- 
tral index with frequency is expected due to energy losses of the 
electrons. 

The frequency scaling of the AME component is the most 
poorly constrained. The distinctive feature is a peak around 20- 
40 GHz (Draine & Lazarian 1998, Dobler & Finkbeiner 2008b, 
Dobler et al. 2009, Hoang et al. 2011). However, a power-law 
behaviour is compatible with most detections above 23 GHz 
(Banday et al. 2003, Davies et al. 2006, Ghosh et al. 2012). This 
could be the result of a superposition of more peaked compo- 
nents along the line of sight or could indicate a peak frequency 
lower than 23 GHz. The most recent WMAP 9-yr results quote 
a peak frequency at low latitudes ranging from 10 to 20 GHz 
for the spectrum in Kr_j units, which means 20-30 GHz when 
considering flux density units. 

3.3. Component separation pipeline 

Several component separation methods adopt the linear mixture 
data model (see Appendix A for full derivation). For each line of 
sight we write: 

x = Hs-\-n, (3) 

where x and n contain the data and the noise signals. They are 
vectors of dimension A^^^ which is the number of frequency chan- 
nels considered. The vector s, having dimension A/c, contains the 
A^c unknown astrophysical components (e.g. CMB, dust emis- 
sion, synchrotron emission, free-free emission, AME) and the 
Nd X A^c matrix H, called the mixing matrix, contains the fre- 
quency scaling of the components for all the frequencies. The 
elements of the mixing matrix are computed by integrating the 
source emission spectra within the instrumental bandpass. When 
working in the pixel domain, Eq. (3) holds under the assumption 
that the instrumental beam is the same for all the frequency chan- 
nels. In the general case, this is achieved by equalizing the reso- 
lution of the data maps to the lowest one. When working in the 
harmonic or Fourier domain, the convolution for the instrumen- 
tal beam is a multiplication and is linearized without assuming a 
common resolution. 

Within the linear model, we can obtain an estimate s of the 
components s through a linear mixture of the data: 

s = Wx, (4) 

where W is called the reconstruction matrix. Suitable reconstruc- 
tion matrices can be obtained from the mixing matrix H. For ex- 
ample: 

\N=[H'^C-'H]-'H^C-' (5) 

is called the generalized least square (GLS) solution and only 
depends on the mixing matrix and on the noise covariance Cn- 

The mixing matrix is the key ingredient of component sep- 
aration. However, as discussed in Sec. 3.2, the frequency spec- 
tra of the components are not known with sufficient precision to 
perform an accurate separation. To overcome this problem, our 
component separation pipeline implements a first step in which 
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Table 2. Summary of ancillary data. 



Central frequency 


Label 


Resolution 


Reference 


[GHz] 




[arcmin] 




0.408 


Haslam 


60 


Haslam et al. (1982) 




Ha 


60 


Dickinson et al. (2003) 






6-60 


Finkbeiner (2003) 


22.8-94 


WMAP 7-yr 


56.8-13.8 


Jarosik et al. (2011) 


94 




60 


Finkbeiner et al. (1999) 


2997 


lOOyum 


5 


Schlegel et al. (1998) 


24983, 2997 


IRIS Band 1,4 


4 


Miville-Deschenes & Lagache (2006) 




E(B - V) 


5 


Schlegel et al. (1998) 



the mixing matrix is estimated from the data and a second one 
in which this result is exploited to reconstruct the amplitudes of 
the components. 

3.3.1 . Estimation of the mixing matrix 

For the mixing matrix estimation we rely on the CCA (Bonaldi 
et al. 2006, Ricciardi et al. 2010), which exploits second-order 
statistics of the data to estimate the frequency scaling of the com- 
ponents on defined regions of the sky (sky patches). We used the 
harmonic-domain version of the CCA, whose basic principles 
of operation are reported in Appendix A. This code works on 
square sky patches using Fourier transforms. It exploits the data 
auto- and cross-spectra to estimate a set of parameters describing 
the frequency scaling of the components. The patch-by-patch es- 
timation prevents the detection of small-scale spatial variations 
of the spectral properties. On the other hand, by using a large 
number of samples we retain more information, which provides 
good constraints, even when the components have similar spec- 
tral behaviour. The CCA has been successfully used to separate 
the synchrotron, free-free and AME components from WMAP 
data in Bonaldi et al. (2007). 

We used a patch size of 20° x 20°, obtained as a trade-off 
between having enough statistics for a robust computation of 
the data cross-spectra and limited spatial variability of the fore- 
ground properties. Given the dimension of the region of interest, 
we have 10 independent sky patches. However, exploiting a re- 
dundant number of patches, widely overlapping with each-other, 
enables us to eradicate the gaps between them and obtain a re- 
sult that is independent of any specific selection of patches. We 
covered the region of interest with patches spaced by 2° in both 
latitude and longitude. By re-projecting the results of the CCA 
on a sphere and averaging the outputs for each line of sight we 
can synthesize smooth, spatially varying maps of the spectral 
parameters (see Ricciardi et al. 2010 for more details). 

3.3.2. Reconstruction of the component amplitudes 

The reconstruction of the amplitudes has been done in pixel 
space at 1° resolution using Eq. (5), exploiting the output of the 
previous step. To equalize the resolution of the data maps, the 
a^m of each map have been multiplied by a window function, 
W^^\ given by a 1° Gaussian beam divided by the instrumen- 
tal beam of the corresponding channel (assumed to be Gaussian 
with full width half maximum (FWHM) as specified in Table 1). 
This corresponds, in real space, to convolution with a beam Bs. 
In order to obtain an estimate of the corresponding noise after 
smoothing, the noise variance maps should be convolved with 
= (Bs)^. We did this again in harmonic- space, after hav- 



ing obtained the window function , corresponding to B^, by 

Legendre transforming squaring the result, and Legendre 
transforming back. 

The smoothing process also correlates noise between differ- 
ent pixels, which means that the RMS per pixel obtained as de- 
tailed above is not a complete description of the noise properties. 
However, the estimation of the full covariance of noise (and its 
propagation through the separation in Eqs. (4) and (5)) is very 
computationally demanding. In this work we take into account 
only the diagonal noise covariance and neglect any correlation 
between noise in different pixels. In a signal-dominated case, 
such as the one considered here, the errors on the noise model 
have very small impact on the results. 

4. AME frequency spectrum 

We modelled the mixing matrix to account for five components: 
CMB; synchrotron emission; thermal dust emission; free-free 
emission; and AME. We neglected the presence of the CO com- 
ponent by excluding from the analysis the 100 and 217 GHz 
Planck channels, which are significantly contaminated by the 
CO lines 7=1^0 and 7 = 2^1 respectively (Planck 
HFI Core Team 2011b). CO is also present at 353 GHz, where 
it can contaminate the dust emission by up to 3% in the region 
of interest, and at 545 and 857 GHz, where the contamination is 
negligible. For the estimation of the mixing matrix we used the 
following dataset: 

- Planck 30, 44, 70, 143 and 353 GHz channels; 

- WMAP 7-yr K band (23 GHz); 

- Haslam 408 MHz map; 

- Predicted free-free emission at 23 GHz based on the Hor 
Dickinson et al. (2003) template corrected for dust absorp- 
tion with the Schlegel et al. (1998) E(B - V) map by assum- 
ing a dust absorption fraction of 0.33. 

We verified that the inclusion of the WMAP Ka-W bands in 
this analysis did not produce appreciable changes in the results. 
The explored frequency range is now covered by Planck data 
with higher angular resolution and sensitivity. Caution is needed 
when using Ha as a free-free tracer: dust absorption (Dickinson 
et al. 2003) and scattering of Hck photons from dust grains (Wood 
& Reynolds 1999, Dong & Draine 2011) cause dust-correlated 
errors in the free-free template, which could bias the AME spec- 
trum. The impact of such biases has been assessed through sim- 
ulations as described in Sect. 4.1. 

For dust emission we used the model of Eq. (1) with = 
18 K and estimated the dust spectral index fid- The reason why 
we fixed the dust temperature is that this parameter is mostly 
constrained by high-frequency data, which we do not include 
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Fig. 3. CCA estimation of the AME frequency spectrum in the region of interest for real data (left panels) and simulated data (right panels). Top: 
areas spanned by the estimated spectra including Icr errors. Middle and bottom: histograms of the spectral parameters m^o and Vp on A^side = 16 
estimated spectral index maps. For the simulated case (right panels) we consider two convex spectra pealdng at 19 GHz and 26 GHz and a power- 
law model. Top right: the true inputs are shown as solid blaclc lines (power-law) with triangles (19 GHz peak) and squares (26 GHz pealc) and 
the estimations as shaded areas. The blue and red colours show estimations done exploiting the free-free templates FFi and FF2 described in 
Appendix C. Middle and bottom right panels: the true inputs are shown by solid and dotted vertical lines for the simulations pealdng at 26 GHz 
and 19 GHz respectively; the blue dot-dashed and red solid histograms show the estimations carried out using the FFi and FF2 templates, and the 
black dashed lines show the estimations for the 19 GHz input spectrum. 



in this analysis. The temperature = 18K is consistent with 
the 1 -component dust model by Finkbeiner et al. (1999) and 
in good agreement with the median temperature of 17.7 K es- 
timated at \b\ > 10° by Planck Collaboration XIX (2011). For 
the dust spectral index we obtained fid = 1.73 ± 0.09. For syn- 
chrotron radiation we adopted a power-law model with fixed 
spectral index JS^ = -2.9 (e.g., IMiville-Deschenes et al. 2008), 
as the weakness of the signal prevented a good estimation of 



this parameter. We verified that diff'erent choices for jS^ (up to a 
10 % variation, fS^ from -2.6 to -3.2) changed the results for the 
other parameters only of about 1 %, due to the weakness of the 
synchrotron component with respect to AME and thermal dust. 
As a spectral model for AIME we adopted the best-fit model of 
Bonaldi et al. (2007), which is a parabola in the log(5' )-log(y) 
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plane parametrized in terms of peak frequency Vp^ and slope at 
60 GHz meo- 

log Trj^ameCv) oc + 2 log y + . . , . . . (6) 

\log(yp/60) / 21og(yp/60) 

Details of the model and justification of this choice are given in 
Appendix B. We also tested a pure power-law model (y oc y^) 
for AME, fitting for the spectral index a, but we could not ob- 
tain valid estimations in this case. This is what we expect when 
the true spectrum presents some curvature, as verified through 
simulations (see Sect. 4.1 and Appendix C). 

Our results for the AME spectrum are shown in the left pan- 
els of Fig. 3. On average, the AME peaks at 25.5 GHz, with a 
standard deviation of 0.6 GHz, which is within estimation errors 
(1.5 GHz). This means we find no significant spatial variations 
of the spectrum of the AME in the region of the sky considered 
here. However, we recall that this only applies to diff'use AME, 
as our pipeline cannot detect small-scale spatial variations, and 
we are restricted to a limited area of the sky. 

Our results on the peak frequency of the AME are similar 
to those of Planck Collaboration XX (2011) for Perseus and 
p Ophiuchi and (Planck Collaboration TBC 2013) on compact 
dust clouds. WMAP 9-yr MEM analysis (Bennett et al. 2012) 
measures the position of the peak for the spectrum in Kr_j units 
and finds a typical value of 14.4 GHz for diff'use AME at low 
latitudes, which roughly corresponds to 27 GHz when the spec- 
trum is in flux density units. According to previous work, at 
higher latitudes the peak frequency is probably lower (see e.g. 
Banday et al. 2003, Davies et al. 2006, Ghosh et al. 2012). 
Interestingly, the same CCA method used in this paper yields Vp 
around 22 GHz when applyied to the North Celestial Pole region 
(towards / = 125°, b = 25°, Bonaldi & Ricciardi 2012). Spatial 
variations of the AME spectrum could explain these differences. 

4.1. Assessment through simulations 

The reliability of our results has been tested with simulations. 
The main purposes of this assessment are: 

- to verify the ability of our procedure to accurately recover 
the AME spectrum for diff'erent input models; 

- to investigate how the use of foreground templates — free- 
free in particular — can bias the results. 

We did this by applying the procedure described in Sect. 4 to 
sets of simulated data, for which the true inputs are known. For 
the first target, we performed three separate simulations includ- 
ing a diff'erent AME model: two spinning dust models, peaking 
at 19 GHz and 26 GHz, and a spatially varying power-law. For 
the second target, we introduced dust-correlated biases in the 
free-free template and quantified their impact on the estimated 
parameters. The full description of the simulations and of the 
tests performed is given in Appendix C. 

The results are displayed in the right panels of Fig. 3. In the 
top panel we show, for each of the three tested input models, the 
true spectrum (solid line) and the estimated spectrum with er- 
rors (shaded area). The red and blue areas distinguish between 
two free-free templates (referred to as FFi and FF2), which are 
biased in a diff'erent way with respect to the simulated free-free 
component. In the middle and bottom panels we show the his- 
tograms of the recovered spectral parameters compared with the 



^ The peak frequency Vp is defined for the specrum in flux density 
units 



true inputs (vertical lines); the red and blue colours are as before. 
We conclude the following: 

- If the input AME is a convex spectrum, we are able to ac- 
curately recover the peak frequency, yp, for both the 19 and 
26 GHz input values. Our pipeline is able to distinguish very 
clearly between the two input models; biases in the free-free 
template do not afl'ect the recovery of the peak frequency. 

- The estimated spectrum can be slightly biased above 40- 
50 GHz, where the AME is faint, as a result of limitations 
of the spectral mode we are using (see Appendix B) and er- 
rors in the free-free template. The systematic error on meo is 
quantified as 0.5-0.6. 

- If the input AME spectrum is a power-law, we obtain a good 
recovery when fitting for a spectral index. 

When the AME is a power-law the parabolic model is clearly 
wrong, as the parameter describing the position of the peak is 
completely unconstrained and the model steepens considerably 
with frequency. Similarly, when the AME is a curved spectrum 
the power-law model is too inaccurate to describe it. As ex- 
pected, both these estimations fail to converge. 

5. Reconstruction of the amplitudes 

The reconstruction of the amplitude of the components has been 
performed on the 1° resolution version of the dataset. We used 
the same frequencies exploited for the estimation of the mix- 
ing matrix, except for the free-free template, which has been ex- 
cluded to avoid possible biases in the reconstruction. The results 
are shown in Fig. 4. The first and second rows show the com- 
ponents reconstructed at 30 GHz (from left to right: synchrotron 
emission, free-free emission, AME, and thermal dust emission) 
and the corresponding noise RMS maps. Thanks to the linearity 
of the problem, the noise variance maps can be obtained by com- 
bining the noise variance maps of the channels at 1° degree res- 
olution with the squared reconstruction matrix W. The noise on 
the synchrotron and thermal dust maps is low compared to that 
for free-free and AME. This is because the 408 MHz map and 
the Planck 353 GHz channel give good constraints on the ampli- 
tudes of synchrotron and thermal dust emission respectively. 

The AME component is correlated at about 60 % and 70 % 
with the lOOyum and the E(B - V) dust templates by Schlegel 
et al. (1998), 40 % with Haslam et al. (1982) 408 MHz and 20 % 
with Ha. This favours emission mechanisms based on dust rather 
than to other hypotheses, such as curved synchrotron emission 
and free-free emission. The E(B - V) template correlates better 
with thermal dust emission than the 100 jim map (the correlation 
coefficients being 0.73 ± 0.01 and 0.96 ± 0.01 respectively). This 
is expected if AME is dust emission. In fact, both spinning dust 
and thermal dust emission are proportional to the column den- 
sity, for which E(B - V) is a better estimator than the lOO^t/m 
emission, which is strongly aff'ected by the dust temperature. 

The errors due to the separation process (third and fourth row 
of Fig. 4) are obtained by propagating via Monte-Carlo the un- 
certainties on the mixing matrix estimated by CCA to the recon- 
struction of the components (see Ricciardi et al. 2010 for more 
details). Essentially, the mixing matrix parameters are random- 
ized according to their posterior distributions; the component 
separation error on the amplitudes is estimated as the variance 
of GLS reconstructions for diff'erent input mixing matrices. 

One complication is that in the present analysis we did 
not estimate the synchrotron spectral index, but we fixed it at 
y^s = -2.9. Thus, we do not have errors on the synchrotron spec- 
tral index from our analysis. We therefore considered two cases: 
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Fig. 4. 1° resolution reconstruction at 30 GHz of (from left to right): synchrotron emission; free-free emission; AME; and thermal dust emission. 
These reconstructions are performed as described in Sect. 3.3.2. Rows from top to bottom: component amplitudes; noise RMS; predicted RMS of 
component separation error due to the estimation of AME and thermal dust spectra; and predicted RMS of component separation error including 
a random error on - -2.9 + 0.1. 



one in which we propagated only the errors on the AME and 
thermal dust spectral parameters, thus assuming no error on 
(third row of Fig. 4); and another in which we included an in- 
dicative random error Ay^s = 0.1 (last row of Fig. 4). 



The predicted error due to separation is generally higher than 
noise and on average of the order 15-20 % of the component am- 
plitude for AME, free-free, and dust. Once we allow some scat- 
ter on the predicted error on synchrotron emission becomes 
of the order of 50 %: this indicates that the reconstruction of this 
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Fig. 5. Left panels: sum of the components vs frequency maps at (from top to bottom) 30, 44, and 70 GHz. The line is the x = j relation. Right 
panels', pixel distribution of the residual (frequency map-sum of the components) maps compared to the best-fit Gaussian distribution. 



component is essentially prior-driven. The inclusion of Ay^s has 
some effect on the error prediction for free-free emission, while 
AME and dust are mostly unaflFected. 

To evaluate the quality of the separation we compared the 
frequency maps with the sum of the reconstructed components 
at the same frequency. In the left panels of Fig. 5 we plot the sum 
of the components for 30, 44, and 70 GHz against the amplitude 
of the frequency map. The comparison is made at 1° resolution 
with A^side = 128 pixels. The dashed line indicates \ht x = y 
relation, which corresponds to the ideal case in which the two 
maps are identical. 



The agreement between data and predictions is in general 
very good. The scatter of the points does not measure the quality 
of the separation but the signal-to-noise of the maps. It increases 
from 30 to 70 GHz, as the foreground signal gets weaker. The 
errors in the component separation show up as systematic depar- 
tures of the data from the prediction. As those are not apparent, 
we also show on the right panels of Fig. 5 the pixel distribution of 
the residual map compared to the best-fit Gaussian distribution. 
At 44 and 70 GHz the scatter, though quite small, dominates the 
residual and covers the systematic effects, with the exception of 
a few outliers, mostly due to compact sources. At 30 GHz the 
scatter is low enough to reveal a feature: a sub-sample of pixels 
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Fig. 6. Estimated spectra of synchrotron emission (dashed line), free- 
free emission (sohd line), thermal dust emission (dash-dotted line) and 
AME (dotted line) for average local properties (top) and for Barnard's 
region (bottom). 



6- Free-free electron temperature 

The intensity of the free-free emission at a given frequency with 
respect to Ha can be expressed as 



rfF(y)[//KR-j] 
HQf[Rayleighs] 



= 14 7 



0.517 . 



^Q0.029/r4 X 1.08 G(y)(y/ 10)-^ (7) 



where G(y) is the Gaunt factor already introduced in Sect. 3.2 
and T\ is the electron temperature in units of 10"^ K. One 
caveat is that dust absorbs Ha emission; this can be quantified by 
the dust absorption fraction /d (Dickinson et al. 2003). Therefore 
/d and Te are degenerate parameters. 

The ratio T^(y)I^Aa can be obtained by comparing the Ha 
and free-free emission from component separation through a 
temperature-temperature plot (TT analysis). We made free- free 
versus Ha plots by using the CCA free-free solution at 30 GHz 
and both Dickinson et al. (2003) and Finkbeiner (2003) Ha tem- 
plates corrected for dust absorption for diff'erent values of /d. 
We considered 3° resolution maps, sampled with A^side - 64 pix- 
els. Besides point sources, we excluded from the analysis the 
region most aff'ected by dust absorption based on the Schlegel 
et al. (1998) E{B - V) map, as shown in the left panel of Fig. 7. 
The electron temperature has been inferred by fitting the data 
points with a linear relation and converting the best-fit slope to 
Te through Eq. (7). The error on has been derived from the 
error on the best-fit slope given by the fitting procedure, through 
error propagation. In the right panel of Fig. 7 we show the TT 
plots for the Dickinson et al. (2003) Ha template, with diff'erent 
colours corresponding to diff'erent values of /d, and the best-fit 
linear relations. The electron temperatures are reported in the 
top part of Table 3. We obtain = 6000-3500 K with U = 0- 
0.5 for the Dickinson template; the Finkbeiner template yields 
higher, but consistent, values (T^ = 7000-4500 K with /d = 0- 
0.5). 



in which the reconstructed signal is higher than the true one, thus 
creating a negative in the residual. 

This kind of systematic eff'ect is very difficult to avoid when 
separating many bright components, because small errors in the 
mixing matrix cause bright features in the residual maps. Our 
Monte-Carlo approach is however able to propagate these er- 
rors. At 30 GHz the brightest components are AME and free-free 
emission, for which the predicted component separation error is 
on average 0.04-0.05 mKcMB, in agreement with the level of the 
non-Gaussian residuals. Coherent structures in the residual maps 
are induced by the low resolution of the maps of spectral parame- 
ters, which means that over nearby pixels the error in the mixing 
matrix, and thus on the separation, is similar. 

In Fig. 6 we show the amplitude of the components as a func- 
tion of frequency. The top panel represents the typical behaviour 
in the Gould Belt, while the bottom one refers to a particular 
case, Barnard's region where free-free emission is particularly 
strong. The points are the average amplitude of the components 
at each frequency within the selected regions of the sky. The 
scaling of the amplitudes with frequency is, by construction, 
given by the spectral model estimated with CCA. The error bars 
measure the scatter induced on the amplitudes by the errors on 
the spectral parameters (also including Ay^s = 0.1). 



6.1. Comparison with cross-correlation witli templates 

An alternative way to compute rff(y)/Ha and Tq is through 
cross-correlation of the Ha template with frequency maps (CC 
analysis). We cross-correlate simultaneously the templates for 
free-free, dust, and synchrotron emission, as described in Ghosh 
et al. (2012). We used the 408 MHz map from Haslam et al. 
(1982) as a tracer of synchrotron emission, Dickinson et al. 
(2003) Ha as a tracer of free-free emission and the Finkbeiner 
et al. (1999) model eight 94 GHz prediction as a tracer of dust 
emission. We used the same resolution, pixel size and sky mask 
adopted for the TT analysis (3° and A^side = 64). As pointed 
out by Ghosh et al. (2012), at this resolution the template- fitting 
analysis is more reliable than at 1° because the smoothing re- 
duces artifacts in the templates. The correlation coefficients are 
computed for each emission process at a given frequency by 
minimizing the generalized;^^ expression. We also fitted for an 
additional monopole term that can account for off'set contribu- 
tions in all templates and the data in a way that does not bias the 
results (Macellari et al. 201 1). The chance correlation of the tem- 
plates with the CMB component in the data causes a systematic 
error in the correlated coefficients and has been estimated using 
simulations. We generated 1000 random realizations of the CMB 
using the best-fit ACDM model ^ and cross-correlated each of 
them using the templates with the same procedure applied to the 
data. The amplitude of the predicted chance correlation, given 
by the RMS over the 1000 realizations, is 1.13/iKcMB//^KcMB 

^ http://lambda.gsfc.nasa.gov/product/map/dr4/pow_tt_spec_get.cfm 
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Table 3. Inferred [K] for TT analysis and CC analysis using different Hof templates and dust absorption fractions /d. 



Method 


Template 


/d = 0.0 


/d = 0.1 


/d = 0.3 


/d = 0.5 


TT analysis 


Dickinson 


5900 + 2500 


5500 + 2500 


4700 + 2200 


3500 + 1600 




Finkbeiner 


6800 + 2500 


6400 + 2500 


5500 + 2400 


4600 + 1900 


CC analysis 


Dickinson 


4900+ 1100 


4000 + 1000 


2600 + 850 


2000+ 700 




Finkbeiner 


7000 + 1300 


6400 + 1200 


5000 + 1000 


3400 + 840 



ap \ \ 1 » 1 0-/- 




-500 ^ . ^ . . ^ . - - 1 

50 100 150 200 

Ha[Rayl9ighs] 

Fig. 7. TT analysis for estimation of T^. Top: gnomonic projection 
showing the mask used (masked pixels are in black, while pixels used 
in the analysis are in white). Bottom: TT plot comparing the CCA free- 
free solution with the Hor template for /d = (points) and linear fits to 
the TT plots for different values of /d (lines). 



for the dust template, 1.12yuKcMB/R for the free-free template 
and 3.8 yuKcMe/K for the synchrotron template. 

In the top panel of Fig. 8 we compare the Yia correlation 
coefficients (points with error bars) with the component separa- 
tion results for free-free emission obtained in Sect. 5 (shaded 
area). The flux for both the component separation and cross- 
correlation has been computed as the standard deviation of the 
maps (the separated free-free map and the scaled Ha template, 
respectively) as this is not affected by possible offsets between 
the Planck data and the Yia template. There is generally good 



agreement between the two results; in the frequency range 40- 
60 GHz there is an excess in the correlated coefficients, which 
could be indicative of a contribution from the AME component 
(similar to that found by Dobler & Finkbeiner 2008b). Flattening 
of the CC coefficients for y > 60 GHz is consistent with positive 
chance correlation between the CMB and the Ha template. 

The dust-correlated coefficients are compared with the com- 
ponent separation results in the bottom panel of Fig. 8. The 
agreement is very good for y < 40 GHz and y > 100 GHz, where 
AME and thermal dust emission are strong. In the 40-70 GHz 
range the CC results are higher than the component separation 
results. As discussed in Appendix B, the parametric fit to the 
AME spectrum implemented by CCA could be inaccurate in 
this frequency range, where the AME is faint. Alternatively, a 
similar effect could be explained by the presence of a secondary 
AME peak, around 40 GHz (e.g. Planck Collaboration XX 201 1, 
Ghosh et al. 2012) or flattening of the dust spectral index towards 
low frequencies, which are not included in our spectral model. 
Discriminating between these hypotheses is not possible given 
the large error bars. 

To determine the free-free electron temperature the Ha 
correlation coefficients have been fitted with a combination 
of power-law free-free radiation (with fixed spectral index of 
-2.14) and a CMB chance correlation term (which is constant 
in thermodynamic units). The amplitude of the free-free com- 
ponent with respect to Ha resulting from the fit, and its un- 
certainty, yield Tq and the corresponding error bar. The results 
for the Gould Belt region outside the adopted sky mask are re- 
ported in the bottom part of Table 3. We find = 5000-2000 K 
for /d = 0-0.5 with the Dickinson et al. (2003) template, and 
Te = 7000-3500 K for /d = 0-0.5 with the Finkbeiner (2003) 
template. 

With respect to the TT analysis, these results are more sen- 
sitive to the choice of the template and the /d correction. In fact, 
the correlation coefficients mostly depend on the brightest fea- 
tures in the map, which can differ significantly from one tem- 
plate to another, while the fit of the TT plot is driven by the 
larger sample of pixels, which is more difficult to modify in a 
systematic way. The CC results are always consistent with the 
TT ones; however, we note that, for the Dickinson template, they 
are systematically lower. One possible explanation is that the 
two methods have different sensitivity to errors in the template 
or in the model used. Spatial variability of Tq could also give 
a systematic difference in the results, since the two approaches 
weight pixels differently. This confirms that estimating the free- 
free electron temperature is a difficult problem and that caution 
is needed when interpreting the results. 



7. AME as spinning dust emission. 

An explanation that is often invoked for the AME is elec- 
tric dipole radiation from small, rapidly spinning. Poly cyclic 
Aromatic Hydrocarbon (PAHs) dust grains (Erickson 1957, 
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Fig. 8. Comparison between correlation coefficients (symbols with er- 
ror bars) and component separation results (shaded areas) for free-free 
emission (top) and dust emission (bottom). For free-free emission we 
show the Dickinson et al. (2003) Ha correlation coefficients and for 
dust emission the Finkbeiner et al. (1999) correlation coefficients. The 
grey area in the bottom panel is the sum of the AME (blue) and thermal 
dust (red) components. The dash-dotted line in both panels shows the 
Icr error due to the chance correlation of CMB with foreground tem- 
plates, estimated using simulations. 



Draine & Lazarian 1998, Dobler & Finkbeiner 2008b, Dobler 
et al. 2009, Hoang et al. 201 1). 

Alternatively, the AME could be due to synchrotron radia- 
tion with a flat (hard) spectral index (e.g. Bennett et al. 2003). 
The presence of such a hard spectrum synchrotron component 
could be highlighted by comparing the 408 MHz map of Haslam 
et al. (1982), which would predominantly trace steep spectrum 
radiation, with the 2.3 GHz map by Jonas et al. (1998), which 
would be more sensitive to flat spectrum radiation. This issue 
has been studied in detail by Peel et al. (2011) using a cross- 
correlation of WMAP 7-yr data with foreground templates. They 
analysed the region defined by 170° < / < 210°, -55° < b < 
-25° and found that the dust-correlated coefficients are mostly 
unaffected by the use of the 2.3 GHz template instead of the 
408 MHz template. This indicates that hard synchrotron radia- 
tion cannot account for most of the dust-correlated component 
at low frequencies. 

To check the hypothesis of spinning dust emission we ap- 
plied the method proposed by Ysard et al. (201 1), which exploits 
the SpDust (Ali-Haimoud et al. 2009, Silsbee et al. 2011) and 
DustEM (Compiegne et al. 2011) codes, to model the frequency 
spectra of thermal and anomalous dust emission from the mi- 



crowaves to the IR. The dust populations and properties are as- 
sumed to be the same as in the diffuse interstellar medium at high 
Galactic latitude (DHGL), defined in Compiegne et al. (2010). 
This model includes three dust populations: PAHs; amorphous 
carbonaceous grains; and amorphous silicates. For PAHs, it as- 
sumes a log-normal size distribution with centroid ao = 0.64 nm 
and width cr = 0.4, with a dust-to-gas mass ratio Mpah/^h = 
7.8 X 10-4. 

By fitting the thermal dust spectrum with DustEM we deter- 
mine the local intensity of the interstellar radiation field. Go (the 
scaling factor with respect to a UV flux of 1 .6 x 10"^ erg s'^cm"^ 
integrated between 6 and 13.6eV), and the hydrogen column 
density, Nu- We then fit the AME spectrum with SpDust, the 
only free parameter being the local hydrogen density nu- We as- 
sume a cosmic-ray ionization rate fcR = 5 x 10"^^ s"^H"\ and 
take the electric dipole moment to be as in Draine & Lazarian 
(1998), a prescription also shown to be compatible with the 
AME extracted from WMAP data (Ysard et al. 2010). It is worth 
noticing that there is a degeneracy with the size of the grains 
(smaller size yields higher peak frequency and intensity of the 
AME). However, the size distribution can only be constrained 
using shorter wavelength data (typically 3-8 yum). The size we 
are adopting (0.64 nm) is motivated by its ability to reproduce 
the data in the mid-IR (Compiegne et al. 2011); other models 
adopt different sizes (e.g. 0.54 nm and 0.5 nm in Draine & Li 
2001 and Draine & Li 2007 respectively). 

As the Gould Belt region contains strong foreground emis- 
sion components, significantly correlated with each-other, we 
expect different environments to be mixed in a complex way. In 
order to obtain meaningful results for the physical modelling we 
tried to isolate sub-regions where single environments dominate. 
To first order, we can use the free-free emission as a tracer of the 
ionized gas environment, CO emission as a tracer of molecular 
gas, and associate the rest of the emission with the diffuse ISM. 
In Fig. 9 we schematically map the different environments by 
setting a threshold on the free-free emission coming from com- 
ponent separation, the CO emission from Planck, and the to- 
tal foreground emission at 30 GHz. We identified two relatively 
big sub-regions (shown as circles in Fig. 9) as selections which 
are dominated by ionized gas and diffuse ISM environments. It 
would not be meaningful to consider smaller areas because of 
the patch-by-patch estimation of the AME frequency scalings, 
which means that our AME spectra are averaged over relatively 
large areas of the sky. Due to the dumpiness of the molecular gas 
environment it was not possible to select a region for this case. 
It is worth noting that some molecular gas may be contained in 
the diffuse ISM region. 

The spectra of AME and thermal dust in the 20-353 GHz 
frequency range are based on component separation results. The 
frequency scaling is that estimated with CCA and the normal- 
ization is given by the average of the reconstructed amplitude 
map in the region of the sky considered. The error bars on the 
data points include the RMS of the amplitude in the same re- 
gion (considered as the error in the normalization) and the er- 
rors on the estimated spectral parameters. The thermal dust spec- 
tra have been complemented with higher frequency data points 
computed directly from the frequency maps: Planck 545 GHz 
and 857 GHz; IRIS 100 yum map; and the IRIS 12yum map cor- 
rected for Zodiacal light emission used in Ysard et al. (2010). 

The results of the modelling for the ionized gas and diffuse 
ISM regions within the Gould Belt are shown in Fig. 10. The 
empirical spectra of AME coming from component separation 
can be successfully modelled as spinning dust emission for both 
regions. The match between data and model becomes worse at 
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Fig. 9. Partitioning of the Gould Belt region based on thresholds over 
free-free emission (red), CO emission (blue), and total emission at 
30 GHz (yellow), used as tracers of H ii gas, molecular gas, and diffuse 
ISM environments, respectively, and the rest of the Gould Belt region 
(light blue). Circled regions are those selected for the computation of 
spectra and have been labelled as Hii-gas (Barnard's arc) and diffuse 
ISM (centred on / = 190°, b - -35°) regions, respectively. 



higher frequencies, where the AME spectrum could be biased 
(see Sect. C and Appendix B). 

The joint fit of thermal and spinning dust models yields 
plausible physical descriptions of the two environments. In the 
top panel of Fig. 10 the diffuse ISM region is modelled with 
A^H = 2.46 X 10^1 Hcm-^ Go = 0.55 and ^H=50cm-l The ion- 
ized region (middle panel) is modelled with A^h = 5.73 x 10^^ 
Hcm-2, Go = 0.90 and nn =25 cm'l 

We tested the stability of these results against calibration er- 
rors on the high frequency Planck (545 and 857 GHz) and IRIS 
(lOOyum and 12yum) data (the remaining data points come from 
the component separation procedure and their error bars already 
include systematic uncertainties). The actual calibration uncer- 
tainty on the Planck 545 and 857 GHz channels is estimated to 
be about 5 %, but is currently being re-evaluated; that on the IRIS 
data is of the order 10% or larger, especially at 12 jim where it 
also includes errors on the zodiacal light subtraction. We have 
verified that very conservative uncertainties up to 20 % both on 
Planck and IRIS data have negligible impact on Go, while they 
may affect A^h and nu (up to a level of about 10 %). The overall 
picture however does not change: the ionized region is less dense 
and illuminated by a stronger radiation field than the diffuse re- 
gion (which is expected to contain mostly neutral gas). Both the 
spectra can be modelled as spinning dust emission arising from 
regions with densities characteristics of the cold neutral medium 
(CNM, a few tens of H per cm^). This confirms the results of 



Planck Collaboration XX (201 1) and Planck Collaboration XXI 
(2011), showing that most of the observed AME could be ex- 
plained by spinning dust in dense gas. In fact, whenever we 
have a mixture of warm neutral medium (WNM), warm ionized 
medium (WIM) and CNM, the spinning dust spectrum is domi- 
nated by the denser phase, which emits more strongly. In the bot- 
tom panel of Fig. 10 we consider for the ionized region a mixture 
of two phases, one having lower density (^h = 0.1 cm"^, 46 %) 
and one having higher density (^h = 55 cm"^, 54%), illumi- 
nated by the same Go as in the middle panel. Such a mixture fits 
the data somewhat better at 23 GHz than the one-phase model 
considered previously (the error in the fit at this frequency being 
0.3 (T instead of 0.9 cr). In order to fully isolate and study dif- 
ferent ISM phases (ionized/neutral, dense/diffuse), both the ob- 
servations and the analysis should be carried out at high angular 
resolution. 



8. Conclusions 

We performed an analysis of the diffuse low-frequency Galactic 
foregrounds as seen by Planck in the Southern part (130° < / < 
230° and -50° < b < -10°) of the Gould Belt system, a local 
star-forming region emitting bright diffuse foreground emission. 
Besides Planck data our analysis includes WMAP 7-yr data and 
foreground ancillary data as specified in Table 2. 

We used the CCA (Bonaldi et al. 2006, Ricciardi et al. 2010) 
component separation method to disentangle the diffuse Galactic 
foregrounds. In the region of interest the synchrotron component 
is smooth and faint, and ultimately not well constrained. 

The free-free emission is strong and it clearly dominates 
in the Orion-Barnard region. We inferred the free-free electron 
temperature both by cross-correlation (CC) of channel maps 
with foreground templates and temperature-temperature (TT) 
plots comparing the reconstructed free-free emissio with Ha 
maps. We obtained Tq ranging from 7000 to 2000 K for /d=0- 
0.5, in agreement with other recent results (Ghosh et al. 2012). 
The use of the Finkbeiner (2003) Ha template yields systemati- 
cally higher Tq than the Dickinson et al. (2003) one. In the case 
of the TT analysis the difference is around 1000 K (< Icr), while 
for the CC analysis it can reach 2500 K (about 2cr). The CC 
results for the Dickinson et al. template are also systematically 
lower than the TT ones, yet consistent within Icr. 

The AME is the dominant foreground emission at low fre- 
quencies over most of the region considered. We estimated the 
AME peak frequency in flux density units to be 25.5 ±1.5 GHz, 
almost uniformly over the region of interest. This is in agree- 
ment with AME spectra measured in compact dust clouds (e.g. 
Planck Collaboration XX 2011) and WMAP 9-yr results at low 
latitudes (once the same convention is adopted, e.g. their AME 
spectrum is converted from Kr_j to flux density, Bennett et al. 
2012). In the case of diffuse AME at higher latitudes a lower 
peak frequency is favoured (Banday et al. 2003, Davies et al. 
2006, Ghosh et al. 2012, Bonaldi & Ricciardi 2012). The abil- 
ity of our pipeline to correctly recover the peak frequency of 
the AME, Vp, has been verified through realistic simulations. We 
also considered the effect of systematic errors in the spectral 
model and in the free-free template and we demonstrated that 
they have negligible impact on Vp. 

Following Peel et al. (2011), a hard (flat spectrum) syn- 
chrotron component would not be sufficient to account for the 
dust-correlated low-frequency emission in this region. In sup- 
port of the spinning dust mechanism, we performed a joint mod- 
elling of vibrational and rotational emission from dust grains as 



13 



Planck Collaboration et al.: Galactic diffuse components in the Gould Belt system 




10^ 10^ 10^ 10* 

Frequency [GHz] 




10^ 10^ 10^ 10* 

Frequency [GHz] 




Frequency [GHz] 

Fig. 10. Frequency spectra (black points with error bars) for thermal 
dust emission and AME compared, respectively, with DustEM and 
SpDust (dashed lines) for the diffuse ISM (top) and ionized gas (middle 
and bottom) regions within the Gould Belt. The solid line is the sum of 
the DustEM and SpDust models. The grey area on the top and middle 
panels correspond to the +1 cr variations on the best- fit riu values, when 
fitting for a single phase. In the bottom panel we consider a mixture 
of two phases (^h = 0.1 cm~^ and riu = 55 cm~^, in the proportion of 
46 % and 54 %, respectively), which marginally improves the fit for the 
ionized gas region at 23 GHz (the error is 0.3 cr instead of 0.9 cr ). 



described by Ysard et al. (2011) and we obtained a good de- 
scription of the data from microwaves to the IR. The fit, which 
we performed separately for the ionized area near to Barnard's 
arc and the diffuse emission towards the centre of our region, 
yields in both cases plausible values for the local density and 
radiation field. This indicates that the spinning dust mechanism 
can reasonably explain the AME in the Gould Belt. 
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Appendix A: Harmonic-domain CCA 

The sky radiation, x, from direction r at frequency v results from 
the superposition of signals coming from A^c different physical 
processes sf 



x(r, v) = ^ Sj(r, y). 

7=1 



(A.1) 



The signal x is observed through a telescope, the beam pattern of 
which can be modelled, at each frequency, as a spatially invariant 
point spread function B(r, v). For each value of v, the telescope 
convolves the physical radiation map with B. The frequency- 
dependent convolved signal is input to an A^d -channel measuring 
instrument, which integrates the signal over frequency for each 
of its channels and adds noise to its outputs. The output of the 
measurement channel at a generic frequency v is 

x^(r) = j B(r- r\ V) ^ ty{V)sj{r\ V)d/dV + ny{r\ (A.2) 

^ 7=1 

where ty{V) is the frequency response of the channel and ^y(r) 
is the noise map. The data model in Eq. (A.2) can be simplified 
by virtue of the following assumptions: 

- Each source signal is a separable function of direction and 
frequency, i.e.. 



Sj{r,v) = Sj(r)fj(y); 



(A3) 



- B(r, y) = By(r) is constant within the bandpass of the mea- 
surement channel. 



These two assumptions lead us to a new data model: 

Xy(r) = By(r) * ^ hyjSj(r) + ny(r), 

7=1 

where * denotes convolution, and 

hyj ^ j ty{V)fj{V)dV . 

For each location, r, we define: 



(A.4) 



(A.5) 



- the A^c-vector s (sources vector) whose elements are Sj{r)\ 

- the A^d- vector x (data vector) whose elements are Xy(r)\ 

- the A^d- vector n (noise vector) whose elements are ^y(r); 

- the diagonal A/d-matrix B whose elements are By{r)\ 

- the A^d X A^c matrix H containing all hyj elements. 



Then, we can rewrite Eq. (A.4) in vector form: 
x{r) = [B * H5](r) + n(r). 



(A.6) 



The matrix H is called the mixing matrix and contains the fre- 
quency scaling of the components for all the data maps involved. 

When working in the pixel domain, under the assumption 
that B does not depend on the frequency, we can simplify 
Eq. (A.6) to 



X = Hs -\- n. 



(A.7) 



where the components in the source vector s are now convolved 
with the instrumental beam. 

Eq. (A.6) can be translated to the harmonic domain, where, 
for each transformed mode, it becomes 



Z = BH5 + A^, 



(A.8) 



15 



Planck Collaboration et al.: Galactic diffuse components in the Gould Belt system 



where X, S, and are the transforms of jc, s, and n, respectively, 
and B is the transform of matrix B. Relying on this data model 
we can derive the following relation between the cross-spectra 
of the data C^(^), sources 0^(0 and noise, C„(0, all depending 
on the multipole ^: 

CA^) = B(^)HC,(OH^B^ (0 + C„(0, (A.9) 

where the dagger superscript denotes the adjoint matrix. 
To reduce the number of unknowns, the mixing matrix is 
parametrized through a parameter vector p (such that H = H (/?)), 
using the fact that its elements are proportional to the spectra of 
astrophysical sources (see Sect. 3.2). 

Since the foreground properties are expected to be spatially 
variable, we work on relatively small square patches of data. 
This allows us to use the 2D Fourier transform to approximate 
the harmonic spectra (see, e.g.. Bond & Efstathiou 1987). 

The HEALPix (Gorski et al. 2005) data on the sphere are 
projected on the plane tangential to the centre of the patch and 
re-gridded with a suitable number of bins in order to correctly 
sample the original resolution. Each pixel in the projected im- 
age is associated with a specific vector normal to the tangential 
plane and it assumes the value of the HEALPix pixel nearest to 
the corresponding position on the sphere. Clearly, the projection 
and re-gridding process will create some distortion in the image 
at small scales and will modify the noise properties. However, 
we verified that this has negligible impact on the spectra in Eq. 
(A.9) for the scales considered in this work and, therefore, on 
the spectral parameters. If x(i, j) contains the data projected on 
the planar grid and X{i, j) is its 2-dimensional discrete Fourier 
transform, the energy of the signal at a certain scale, which cor- 
responds to the power spectrum, can be obtained as the average 
of Z(/, j)X\i, j) over annular bins D^, i = 1, • • • , ^max (Bedini & 
Salerno 2007): 

=TfYj X(iJ)X\iJ), (A. 10) 

where is the number of pairs (/, j) contained in the spectral 
bin denoted by D^. Every spectral bin £ is related to a specific £ 
in the spherical harmonic domain by 

^ = (^-l)2M^/A^pix (A.ll) 

where p is the thickness of the annular bin, = 180/Ldeg (A^pix- 
1), and Ldeg, A/pix are the size in degrees and the number of pixels 
on the side of the square patch, respectively. 

If we reorder the matrices Cx(i) - C„(l) and Cs(^) into vec- 
tors d(£) and c(£), respectively, we can rewrite Eq. (A.9) as 

d(i) = HMc(£)^eih (A. 12) 

where Hk(i) = [B(^)H] [B(^)H], and the symbol denotes the 
Kronecker product. The vector d(f) is now computed using the 
approximated data cross-spectrum matrix in Eq. (A. 10) and €(£) 
represents the error on the noise power spectrum. 

The parameter vector p and the source cross-spectra are fi- 
nally obtained by minimizing the functional: 

<i>[p,cv]= (A.13) 
[dv - HkBip) • Cvf^^B^^v - ^kBip) • cv] + AclCcv. 

The vectors dy and Cy contain the elements d(i) and c(i), re- 
spectively, and the diagonal matrices HkB and the elements 
Hk(i) and the co variance of error 6(i) for all the relevant spectral 



bins. The term ACyCcy is a quadratic stabilizer for the source 
power cross-spectra: the matrix C is in our case the identity ma- 
trix, and the parameter A must be tuned to balance the eff'ects of 
data fit and regularization in the final solution. The functional 
in Eq. (A.13) can be considered as a negative joint log-posterior 
for p and Cy, where the first quadratic form represents the log- 
likelihood, and the regularization term can be viewed as a log- 
prior density for the source power cross-spectra. 

Appendix B: Spectral model for AME 

Theoretical spinning dust models predict a variety of spectra, 
which can be substantially diff'erent in shape, depending on 
a large number of parameters describing the physics of the 
medium. The number of such physical parameters is too large 
to be constrained by the data in the available frequency range. 
For the purpose of the estimation of the spectral behaviour of the 
AME we adopt a simple formula depending on only a few pa- 
rameters. The CCA component separation method used in this 
work implements the parametric relation proposed by Bonaldi 
et al. (2007) [Eq. (6)], depending on the peak frequency, Vp, and 
slope at 60 GHz, m^o. To verify the adequacy of this parametriza- 
tion we produced spinning dust spectra for diff'erent input phys- 
ical parameters with the SpDust code and fitted each of them 
with the proposed relation by minimizing the for the set of 
frequencies used in this work. The input models we consider are: 
weak neutral medium (WNM); cold neutral medium (CNM); 
weak ionized medium (WIM); and molecular cloud (MC). Both 
the input SpDust parameters and the best-fit meo, Vp parameters 
for each model are reported in Table B.l. For comparison, we 
also consider alternative parametric relations and in particular: 

- the model implemented in the Commander component sepa- 
ration method (Pietrobon et al. 2011, Planck Collaboration 
IX 2012) which is a Gaussian in the Tqmb - ln(v) plane, 
parametrized in terms of central frequency and width; 

- the Tegmark et al. (2000) model, which is a modified black- 
body relation [Eq. (1)] having temperature around 0.25 K 
and emissivity index around 2.4; 

As this test does not account for the presence of the other com- 
ponents and does not include any data simulation, it verifies the 
intrinsic ability of the parametric model to reproduce the actual 
spectra. Realistic estimation errors for the CCA model are de- 
rived through simulations in Appendix C. 

Fig. B.l compares the input spectra with the best-fit models 
for the diff'erent parametrizations. In general, the fits are accu- 
rate at least up to v = 50-60 GHz, while at higher frequencies 
the parametric relations may not be able to reproduce the input 
spectra in detail. This is a consequence of fitting complex spec- 
tra with only a few parameters. The fit tends to fail where the 
AME signal is weaker. 

Over the frequency range considered, CCA and Commander 
models fit the input spectrum generally better than the Tegmark 
et al. (2000) model (which falls off' too rapidly at high frequen- 
cies). When adding lower frequency data, however, CCA and 
Commander models will be increasingly inaccurate, as they are 
symmetric with respect to the peak of the emission. The mod- 
els implemented by CCA and Commander perform quite simi- 
larly, despite the diff'erent formulation. As a result, those meth- 
ods are able to give consistent answers, which ensures con- 
sistency between different analyses within Planck (e.g. Planck 
Collaboration IX 2012). 

The CCA model used in this work provides a reasonable 
fit to theoretical spinning dust models for a variety of physical 
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Fig.B.l. Theoretical spinning dust models produced with SpDust 
(solid lines) and fitted with CCA (triangles), Commander (diamonds), 
and Tegmark et al. (2000) (asterisks) models. Input SpDust parameters 
and best-fit parameters for the CCA model are provided in Table B.l. 



conditions. The best-fit parameters that we obtain, reported in 
Table B.l, vary significantly from one input model to another 
and have a straightforward interpretation in terms of the spec- 
trum. 



Appendix C: Description of the simulations 

We simulated Planck and WMAP 7-yr data by assuming 
monochromatic bandpasses positioned at the central frequency 
of the bands, Gaussian beams at the nominal values indicated in 
Tables 1 and 2, and Gaussian noise generated according to real- 
istic, spatially varying noise RMS. Our model of the sky consists 
of the following components: 

- CMB emission given by the best-fit power spectrum model 
from WMAP 7-yr analyses; 

- synchrotron emission given by the Haslam et al. (1982) tem- 
plate scaled in frequency with a power-law model with a spa- 
tially varying synchrotron spectral index as modelled by 
Giardino et al. (2002); 

- free-free emission given by the Dickinson et al. (2003) Ha 
corrected for dust absorption with the E{B - V) map from 
Schlegel et al. (1998) with a dust absorption fraction /d = 
0.33, and scaled in frequency according to Eq. (2) with = 
7000 K; 

- thermal dust emission modelled with the lOOyum map from 
Schlegel et al. (1998), scaled in frequency according to 
Eq. (1) with Td = 18 K and a spatially varying y^d having 
average value of 1.7; 

- AME modelled by the E(B - V) map from Schlegel et al. 
(1998) with intensity at 23 GHz calibrated using the results 
of Ghosh et al. (2012) for the same region of the sky. 

We adopted more than one spectral model for the AME. We 
first considered two convex spectra, generated with the SpDust 
code: one peaking around 26 GHz and the other peaking around 
19 GHz. We also tested a spatially varying power-law model 
(with spectral index of -3.6 ± 0.6, Ghosh et al. 2012), which 
could result from the superposition of multiple convex compo- 
nents along the line of sight. 

It is worth noting that the simulated sky is more complex 
than the model assumed in the component separation. This has 



been done intentionally, to reflect a more realistic situation. 
Another realistic feature we included is the presence of errors 
in the synchrotron and free-free templates. The spatial variabil- 
ity of the synchrotron spectral index modifies the morphology 
of the component with respect to that traced by the 408 MHz 
map from Haslam et al. (1982). The use of Ha as a tracer of 
free-free emission is aflTected by even larger uncertainties. Our 
uncertainties on the dust absorption fraction /d (estimated to be 
/d = 0.33^q|^ at intermediate latitudes by Dickinson et al. 2003) 
and on the scattering of Ha photons from dust grains, can cre- 
ate dust-correlated biases in the template. This is illustrated in 
Fig. C.l, where we compare diflferent versions of the free-free 
template. FFref is our reference template, adopted for the anal- 
ysis of real data and for simulating the component, which is cor- 
rected for /d = 0.33 as described in Dickinson et al. (2003). Two 
more templates (FFi and FF2) have been obtained by correcting 
Ua for /d = 0.33 - 0.15 and /d = 0.33 + 0.1 (±lcr according 
to Dickinson et al. 2003). A final template (FF3) has been ob- 
tained by correcting FFref for scattered light at the 15 % level 
by subtracting from the free-free map the 1 yum map of Schlegel 
et al. (1998) multiplied by a suitable constant factor (Witt et al. 
2010). DiflTerence maps (FFi - FFref) /FFref, presented in the 
lower panels of Fig. C.l, are on average of order of 10%, but 
can be much higher (up to 50-60 %) in regions of strong dust 
emission. 

When analysing the simulated data, we used both FFi and 
FF2 as free-free templates in place of FFref, which corresponds 
to the simulated component. For synchrotron emission the mor- 
phological mismatch between the simulated component and the 
template has been achieved by scaling the component from 
23 GHz to 408 MHz with a spatially varying spectral index. The 
comparison between component and template is presented in 
Fig. C.2; the diflTerences are of the order of 10%. The sim- 
ulated data-sets described above have been analysed with the 
CCA method using the same procedure applied to the real data; 
the results of this assessment are presented in Sect. 4.1. As a 
separate test, we verified the impact of the CMB component on 
the results for Vp and m^o. We generated 100 sets of mock data 
having the same foreground emission and diflTerent realizations 
of CMB and instrumental noise, and repeated the estimation of 
the AME frequency scaling. For this test we used the simulation 
with a spatially constant AME spectrum peaking at 26 GHz. As 
this analysis is computationally demanding, the CCA estimation 
has been performed only on the 10 independent patches covering 
the Gould Belt region (centred on latitudes -20° and -40° and 
longitudes of 140°, 160°, 180°, 200°, and 220°). In Fig. C.3 we 
show the average (diamonds) and RMS (error bars) Vp and m^o 
over the 100 realizations for each patch, for different patches on 
the X-axis. The scatter between the results obtained for diflTerent 
patches (indicated by the grey area in the plots) is typically larger 
than the error bars, measuring the scatter due to different CMB 
realizations. This means that the foreground emission generally 
dominates over the CMB as a source of error. Larger error bars 
associated with the CMB are obtained for three patches having 
fainter foreground emission. For such patches the estimated er- 
rors on Vp and m^o are consistently larger. The CMB variation re- 
sults on average in Avp = 0.1 GHz and Ameo = 0.3, which reach 
0.3 GHz and 0.8, respectively, for the worst sky patch. Those 
values are below the error bars resulting from the analysis of the 
data, which amount to 1-1.5 GHz for Vp and 1.5-2 for meo. The 
CMB has limited impact on the results because this component is 
modelled in the mixing matrix. Having a known frequency scal- 
ing, the statistical constraint used by CCA is able to trace the 
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Table B.l. Spdust input parameters for the spectra in Fig. B.l and best- fit parameters for the CCA spectral model. 





SpDust 










CCA 




model name 


nH [cm-3] 


T[K] 


X 






Vp 


^60 


WNM 


0.4 


6000 


1.00 


0.10 


0.0003 


24.22 


7.53 


CNM 


30.0 


100 


1.00 


0.0012 


0.0003 


29.00 


4.93 


WIM 


0.1 


8000 


1.00 


0.99 


0.001 


27.30 


5.66 


MC 


300 


20 


0.01 


0.0 


0.0001 


38.77 


2.00 






Fig. C.l. Free-free templates at 23 GHz used for the analysis. The reference template FFref is in the upper left comer; the other columns (left to 
right) are FFi, FF2 and FF3, respectively. The diff'erences in the lower panels (FFi - FFref)/FFref are on average of the order of 10 %, but reach 
50 % in regions of strong dust emission. 



pattern of the CMB through the frequencies with good precision 
and hence identify it correctly. 
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Fig. C.2. Upper panel: simulated synchrotron component (left) and syn- 
chrotron template (right) at 23 GHz. Lower panel: difference map di- 
vided by the simulated component. 
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Fig. C.3. Average and RMS of Vp and m^o estimated over simulations 
having different CMB and noise realizations for different patches on 
the X-axis. The grey area is the average and RMS over different patches, 
which is typically larger than that due to noise and CMB. 
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